Particle renormalizations in presence of dissipative environments 
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We study the Aharonov-Bohm oscillations of a charged particle on a ring of radius R coupled to a 
dirty metal environment. With Monte-Carlo methods we evaluate the curvature of these oscillations 
which has the form 1/M* R 2 , where M* is an effective mass. We find that at low temperatures T the 
curvature approaches at large R > I an R independent M* > M, where I is the mean free path in 
the metal. This behavior is also consistent with perturbation theory in the particle - metal coupling 
parameter. At finite temperature T we identify dephasing lengths that scale as T -1 at R > I and 
as T~ 1/4 at R < I. 

PACS numbers: 73.43.Nq,73.23.Ra,74.40.+k 
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I. INTRODUCTION 

The problem of interference in presence of a dissipative 
environment is fundamental for a variety of experimental 
systems. Interference has been monitored by Aharonov- 
Bohm (AB) oscillations in mesoscopic rings 1,2,3 or in 
quantum Hall edge states^ in presence of noise from gates 
or other metal surfaces . Cold atoms trapped by an atom 
chip are sensitive to the noise produced by the chip£££ 
In particular giant Rydberg atoms are studied^ whose 
huge electric dipole is highly susceptible to such noise. 

An efficient tool for monitoring the effect of the en- 
vironment, as proposed by Guinea 9 , is to find the AB 
oscillation amplitude as function of the radius R of the 
ring. This amplitude is measured by the curvatur e) 10 ' 11 ? 12 
of the ground state energy Eq at external flux <j) x = 0, 
i.e. 1/M*R 2 = d 2 E /d(jji1\o, defining an effective mass 
M*. For free particles of mass M this curvature is the 
mean level spacing 1/MR 2 . The particle can be coupled 
to a variety of environments, with three systems of par- 
ticular interest: (i) a Caledeira-Legget (CL) bath 9 , (ii) 
a charged particle in a dirty metal environment*^ 3 - and 

(iii) a particle with an electric dipole in a dirty metal 
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A non-equilibrium study^ has found dephasing lengths 
that have the same power laws, establishing a connection 
between equilibrium and non-equilibrium results. 



II. THE MODEL 

The time dependent angular position 6 m (r) of a par- 
ticle on the ring has in general a winding number m so 
that m (r) = 9(t) + 2tttoTt where 0(0) = 0(1/T) has 
periodic boundary condition. In presence of an external 
flux (f> x (in units of the flux quantum hc/e) the partition 
sum has the form 



Z = e 2lxm ^ J V6e 



= WR 2 / I tt- + 2irmT ) dr 



o \9t / 

l/T ,1/T J T 2] 



JO 
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environment^. System (i) has been studied with a large where the effect of environments, in each of the 3 cases, 
variety of methods, all showing that the AB amplitude is 
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exponentially suppressed ~ e T lR , i.e. a new length 
scale ~ 1/Vr is generated by the coupling 7 to the K( z ) 
environment 9 . System (ii) has been studied by renor- 
malization group (RG) methods*^ finding M* ~ R^ 
with a small /Lt, a Monte Carlo (MC) numerical method 
gavei 3 - /1 — 1.8 at sufficiently large R, while a variational 
scheme^ gave (x = 0. System (iii) was also studied within 
the variational scheme^, leading to /x = as well. 

In the present work we use MC methods to analyze 
mostly system (ii). We find that the energy cutoff used 
in a previous studyi 3 - is insufficient and a higher cutoff 
lo c is needed. In particular we find that at large R > / 
the effective mass M* is R independent, i.e. /z = 0, 
where £ is the mean free path in the metal. For R > i 
we also find that at temperature T the data scales as 
TR, identifying a length scale — l/T. For R < I the 
system reduces to a CL one with a ~ J 1 — 1/4 length scale. 
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Case (i) is the CL system where 7 is the coupling to a 
harmonic oscillator bath; case (ii) is a charge coupled to 
a dirty metal where kp is the Fermi wavevector, I is the 
mean free path in the metal, and r — R/l; case (iii) is an 
electric dipole of strength p coupled to a dirty metal. 

We note that the forms (ii) and (iii) are basedi 3 ^ 4 on a 
wavevector and frequency dependent dielectric function 
for the metal of the form e(q, uj) — 1 + 4ira/ (—iuj + Dq 2 ) 
valid at q < l/£, where a is the conductivity and D is the 
diffusion constant of the metal. The q integrals are cutoff 
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by q < 1 /£ , hence the the forms (ii) and (iii) are valid at 
r > 1. We will use below these forms also at r < 1 since 
they represent qualitatively the decrease of K{z) with r. 
Furthermore, at r — > the form (ii) reduces to that of 
the CL model (i) with aci = 2ar 2 . 

We also note that in model (ii) a < 1 for relevant 
metals. However, model (iii) allows for a large a since the 
dipole parameter p can be large, as e.g. in the Rydberg 
atoms^. 

We are interested in the effect of the environment on 
the visibility of quantum interference as measured by the 
particle. As a measure of this visibility we consider the 
curvature of the Aharonov-Bohm oscillations 
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where F = — T In Z. It is useful to consider a free particle 
a = 0. for which 
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M*(T)J a=Q 
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where t = 2MR 2 T. This identifies the thermal length 
L T ~ 1/VMT. 

In the interacting system a high energy cutoff can be 
identified by considering r — * r' (corresponding to high 
frequencies oj) so that expansion of K(z) and Fourier 
transform yield 



S (r. 



^[MR 2 u? + 2naK" '(O)M]|0(w)| s 
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(27rm) 2 [iMi? 2 T + aK"(0)] . 
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The term linear in |cj| is typical for dissipative systems, 
i.e. the environment induces dissipation on the particle. 
The cutoff uj c is now identified when the kinetic ~ uj 2 
and \lo\ interaction terms are comparable, i.e. 



2iraK"(0) 
MR 2 
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This oj c replaces a possibly higher environment cut- 
off, since significant renormalizations start only below 
lo c where the linear \w\ dispersion leads to In a; terms 
in perturbation theory and to the need for either RG 
treatment, or an equivalent variational scheme 1 ^. Note 
that K"(0) = ^;r 2 ;3r 2 in the 3 models above, hence 
= 7T7/M in case (i), while ui c ~ a /Ml 2 in cases (ii) 
and (iii). 



III. MONTE CARLO PROCEDURE 

For the MC numerical method we need to discretize 
the time axis into a Trotter number Nt of segments, i.e. 
the time interval of each segment is At = 1/(TNt). The 
discrete action is 



UMR 2 N T T + aK"(0)] 
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sin (7r(n - n')/N T ) 



(7) 



The ^aK"(0) term comes from the n = n' interaction 
term by expanding K(z) around z = 0. A key issue in 
our MC study is the choice of energy cutoff 1/Ar and 
the corresponding Trotter number Nt = l/(TAr). The 
correct choice is such that the free kinetic term dominates 
over the single n = n' interaction term, i.e. Nt ^ w c /T, 
with lu c from Eq. ([5]). Hence At rs l/u> c corresponds to 
the cutoff uj c as identified by RG or variational methods. 
A previous MC study on the charge problem^ has chosen 
Nt in the range 1/t to 4/i, i.e. an energy cutoff of w 
l/MR 2 . For large r this cutoff is much smaller than oj c 
and is therefore insufficient. 

Eqs. (dEl identify 1/M*(T)R 2 = 2n 2 T{m 2 )\ (j)x=0 so 
that the MC evaluates the fluctuations in winding num- 
ber (to 2 ) at external flux <j> x = 0. The procedure is to 
start with some to, update 9 n at a time position n to 9 n 
and accept or reject the change according to the MC rule 
with probability exp[S im ^{9 n }-S^{e' n }}. After the N T 
points are successively updated, the winding number is 
shifted to m' = m ± 1 and the shift is accepted or re- 
jected with the probability exp[S'( m ){6' n } - S^ m '^>{9 n }}. 
An update of 9 n is done randomly with a step size that 
produces an acceptance ratio of about 50%ii. 

The inset in Fig. 1 shows the Nt dependence of M/M* 
for the charge problem with r = 5, t = 0.2, a — 0.019. A 
choice for Nt in the range — 4/i is clearly insufficient; 
saturation sets in around Nt ~ 100 which is of order of 
uj c /T — 30. In the following we choose our Nt, in the 
charge problem, to be Nt = 40or 2 jt — 10u) c /(irT), i.e. 
Nt = 95 for the inset parameters. For the dipole case, 
where ui c is 3 times higher we choose Nt = 120ar 2 jt — 
10cj c /(ttT). Fig. 1 shows that for r = 5, t = 0.2, a = 0.02 
(red squares) saturation indeed sets in near Nt — 300. 

This high value of Nt restricts realistic MC studies. 
We have noticed, however, that this high Nt is necessary 
only in the vicinity of n = n' in the double sum of l[7]). 
where the summand is rapidly varying. Hence the double 
sum is taken over all points only in the vicinity of the 
singularity, i.e. for \n — n'\ < 0.03Nt- For points that 
are further separated we coarse grain the sum with fewer 
points, corresponding to an effective Nt = 1/t. 

The results of this procedure are shown by the green 
circles in Fig. 1, and are in agreement with the full calcu- 
lation that includes all Nt points. The double sum has 
then ps il0~ 3 A^ + \t~ 2 terms, much less then the \N T 



terms of the full calculation. We also show data where 
the double sum is coarse grained at all points, including 
those near n = n', by blue triangles. Here the double sum 
has only \t~ 2 terms; this data has significant deviations 
from the full calculation. 

We proceed to discuss our error estimates. At low tem- 
peratures we evaluate (to 2 ), and the average involves typ- 
ically many values of to. To estimate errors we evaluate 
the correlation function for a given run and deduce a 
correlation length £. We discard the initial 10 4 MC iter- 
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FIG. 1: Trotter number dependence of the effective mass for 
the dipole case with r = 5, t = 0.2, a — 0.02, using (i) all Nt 
points in the double sum Eq. (O - red squares, (ii) For points 
\n — n'\ > O.OSNt sum is coarse grained (see text) - green 
circles, (iii) the whole sum is coarse grained - blue triangles. 
Inset: The charge case with r = 5, t = 0.2, a = 0.019 using 
all Nt points in the sums. 



ations and then evaluate the standard deviation a of the 
average data; the error is thenii oV2£ + 1. We typically 
find a short correlation length of a few units and we run 
till an error of ~ 2% is achieved; the number of iterations 
is then w (1 — 2) • 10 5 and in some cases up to 10 6 , where 
each iteration is an update of Nt values of the 9 n . 

At high temperatures t > 1, where M/M* < 10 ~ 3 , 
the probability of to 7^ becomes extremely small so 
that just to = ±1 determine the outcome^. Hence we 
evaluate (m 2 ) = 2(e s ' 1 ~' So - ) )o, averaging with e~ s ". In 
this method we find a rather long correlation length of 
~ 10 3 , yet there is no need to vary m and a 2% accuracy 
can be achieved after w (1 — 2) • 10 5 iterations. 



FIG. 2: AB curvature as function of reduced temperature 
with a = 0.019. All r > 3 vlaues fit the renormalized form 
0.9/(t/0.9) - the lower curve. At r < 1 the data approaches 
f(t) of a free particles - the upper curve. 
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IV. MC RESULTS 

We present here our data for the dirty metal, system 
(ii). In Fig. 2 we show our data for a — 0.019 at low 
temperatures, t < 0.3; we note saturation at t < 0.2. In 
Fig. 3 we collect the limiting low t values of our data for 
various alpha, typically achieved at t w 0.1 — 0.01. The 
data is limited to Trotter numbers N T = 40ar 2 ft < 9000. 

We compare in Fig. 3 the data with results of pertur- 
bation theory (Appendix I). The perturbation is formally 
first order in a, however, it should be valid also for large 
a and small r such that x < 2, where at t = we define 
x = M*{t = 0)/M. The perturbation curves are a good 
fit to the data for r < 1, while at r > 1 and small a the 



FIG. 3: t = limiting values of x = M*(t = 0)/M for var- 
ious a. The full lines are results of perturbation expansion 
(Appendix I). 



fit is qualitatively good, in the sense that saturation is 
achieved at large r. We have also attempted to fit these 
data by a scaling function of the form x = \ + r 2 ~ c g(ar°), 
that is consistent with the r — > form of the perturba- 
tion expansion. In particular, this form with c = 2 would 
scale onto the CL system at r — ► 0, a — ► 00. However, we 
could not find a satisfactory fit even for the small r < 2 
regime. 

Our data shows for the lowest a — 0.019 and for r > 3 
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FIG. 4: AB curvature including high temperatures with a = 
0.019. All data fall in between the upper line f(t) and the 
lower line 0.9/(t/0.9). 

that M/M* reaches saturation with M/M* w 0.9, al- 
most independent of r. The data at r = 20 (shown in 
Fig. 2) is consistent with this saturation, though it is 
not shown in Fig. 3 to keep a convenient scale. In view 
of this saturation at 3 < r < 20 we expect it to persist 
at higher r. In terms of M* ~ r M , our data shows that 
fi < 0.05 and is consistent with // = 0. We note that 
with our revised values of Nt we were not able to reach 
a saturation regime at larger a, see Fig. 3. 

Our result shows that the AB curvature ~ 1/R 2 is 
the same as for free particles, i.e. the ground state has 
no anomaly, at least for weak a — 0.019. Furthermore, 
Fig. 2 shows that M* determines the finite tempera- 
ture behavior, as long as T « w c . Thus if we replace 
M — > M* = M /0.9 in Eq. (gj we obtain the lower curve 
0.9/(^/0.9) in Fig. 2 which is a good fit to the data, The 
thermal length is then L T ~ 1/VM*T. 

In Fig. 4 we show our r > 3 data up to t — 2. The 
data falls in between two lines: 0.9/ (t/0.9) and f(t). The 
lower curve 0.9/(i/0.9) corresponds to the renormalized 
system and fits data with T <C co C) i.e. t <C 4irar 2 . For 
a fixed t as r decreases T approaches lo c and the data 
approaches the upper curve which is the unrenormalized 
free particle form f(t). 

We therefore parameterize our data by a function 
x(r,t) such that M/M* = f(tx)/x. In this way we avoid 
the obvious t dependence associated with mass renor- 
malization and focus on additional temperature effects. 
In Fig. 5 we show that for r > 1 the data for x(t, r) 
scales with t/r. Since t ~ TR 2 the scaling parame- 
ter is ~ TR, identifying a length scale ~ 1/T. A de- 
phasing length scale has been recently derived in a non- 
equilibrium studyi^ which for r > 1 indeed scales with 
1/T. We propose therefore that the additional T de- 



FIG. 5: Scaling of the x variable in M/M* = f(tx)/x for 
r > 1 cases, with a = 0.019 and a — 0.057. 

pendence embedded in our variable x(t, r) is related to 
dephasing of the non-equilibrium situation. 

We note that the perturbation expansion yields for r 

1, 

5fT = l-4«+0(f l*r) r»l. (8) 

While the dependence on t/r is consistent with Fig. 5 
(up to a In r factor), we note that the t/r form in the 
perturbation form ([5]) is valid only at t <C 1 and r > 10. 
Hence the observed scaling, Fig. 5, with t/r up to t « 1 
and at 3 < r < 20 is an unexpected feature. 

In Fig. 6 we show that for r«l the data scales as tr 2 . 
At tr 2 < 0.04 both x(t, r) and a;(0, r) are close to 1 and 
the errors in l/x(t, r) — l/a;(0, r) are too large to draw a 
conclusion in this regime. The same difficulty is with all 
data of small a, hence Fig. 6 shows only a — 0.2, 1. At 
tr 2 > 0.04 the data in Fig. 6 supports a tr 2 scaling. Since 
t ~ TR 2 this implies a length scale ~ T -1 / 4 . We note 
again that similar dependence for a dephasing length was 
found for r < 1 in the non-equilibrium studyi^. 

For r <C 1 we can use the perturbation result Eq. 
(EH 

— - l-2a^ a „ + 4tar 2 r«l. (9) 

n 

This shows the ar 2 scaling at tar 2 < 1. It is remarkable 
that our data in Fig. 6 supports ar 2 scaling up to rather 
high temperatures of t < 1. 

As noted above, the r dependence of K (z) is reliable 
only at r > 1 where the low q, u> form of e(q, to) can be 
used, or at r < 1, which is the CL limit. In fact, for 
a general e(<?, u>) one can expand the response in R and 
obtain that the leading term is K (z) ~ R 2 , i.e. the CL 
form. We conclude then that at both small and large r, 
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ground state&i 3 - of metals. Previous studies proposed 
M* ~ r M with either^^ a small /i or 13 /x = 1.8 or— 
li = 0. Instanton based arguments suggested^ 3 - a M*(r) 
dependence for ar > 1. 

With our revised values of iVy we were able to reach 
a reasonably large r only for weak coupling, a — 0.019. 
For this coupling we observe saturation at 3 < r < 20. 
Although we cannot strictly rule out /j ^ at higher 
r, we find it highly unlikely that an r dependence will 
reappear at r > 20. We propose then \i = at a — 0.019, 
implying /i = at all a (if larger a would show a /i =/= 
it would imply an ulikely singular line in the a, r plane). 
We propose then that — for all a at r 3> 1 and that 
the effect of the environment is a mass renormalization, 
in agrement with the variational studyi^. 



FIG. 6: Scaling of the variable ^| r ^ — x ^ ^ for r <C 1 cases, 
with a = 0.2 and a = 1. 

where if (z) is reliable, the T dependent length scale of 
the equilibrium observable M* /M can be identified with 
a dephasing length. 

V. DISCUSSION 

The possible dependence of M*(r) at T = has been 
of interest as a means of monitoring anomalies in the 



We have found temperature dependent length scales. 
For r > 1 we find T" 1 , while for r < 1 we find T -1 / 4 . We 
note that the same T dependence was found for dephas- 
ing lengths in a nonequilibrium study based on the purity 
of a reduced density matrix^ for the dirty metal situa- 
tion. A dephasing length was deduced^ by comparing a 
dephasing rate with a mean level separation as a condi- 
tion for coherence. It is remarkable that the agreement in 
these dephasing lengths is obtained in both regimes r > 1 
and r <C 1 where the form of Eq. is valid for a dirty 

metal environment; the r <^ 1 form is also valid for other 
realizations of a CL environment. We have therefore the 
intriguing observation that equilibrium scales can iden- 
tify non-equilibrium dephasing length scales. 



APPENDIX A: PERTURBATION EXPANSION 



Consider the action of a particle on a ring in presence of a dissipative environment and a flux <j) x through the ring 
Eq. ([T]) with the dirty metal environment: 

K(z) = 1 - [4r 2 sin 2 | + ip 1/2 = £ a n sin 2 (inz); a = (Al) 

71—1 F 

For a low T expansion it is efficient to perform a duality transformation using the Poisson sum: 

/oo 
#£e 2 ^ 5 (0) (A2) 

m -°° p 

where the sums m,p run on all integers. Hence Eq. ([1]) becomes 

/oo 
d^^e 27 ^^^ 2 * 02 x [1 - 

a £ a n f dr f dr' , 2 - (1 - cos(2™T0(t - r')) (cos[n(fl(r) - 9(r'))}) )} 

^ Jo Jo 2 sin [ttT{t - t')\ 

(A3) 
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where t = 2MR 2 T, j3 = 1/T, Z\ = J T>9 exp(— Si{6}) and the (...) average is taken with respect to exp(— Si), where 

Si{6} = J\t\MR 2 (^ (A4) 

For a Gaussian average we have 

(cos[n(0(T) - 0(t'))]) o = exp[-In 2 ((d(r) - 9(r')) 2 ) } 

= exp[-^ £ 1 ' C0S ^ T ' T/) ] = e-V-r'l/^ (A5) 

where 0(t) = ^ ^ w e _lwr 6'(a;) and w are Matsubara frequencies u> = 2irTx integer. 

For periodic functions we can change integration variables to t\ = r — r', t 2 = ~(r + r') with J" dr2 = /?, and |ti| 
in (|A5[) is chosen in the range (— /3/2,/3/2) to allow for periodicity and continuity at t\ = 0; hence, 



2n-i0(p+^)-7r 2 t0 2 



Z = Zi / d(j) ' 

J-oo p 

/•/3/2 2 T 2 

[l-Pa^n (l-cos(2^r0r)e-" M/*)] (A6) 

^ 7^/3/2 2sm (vrTr) 

Integrating we obtain 



E. _l£±*xli „ \ - f 13 / 2 TT 2 T 2 , (p + .»x) 2 . (p + .» x - n T T ) 2 „-|: 



. iP+fet-Tr)- " 2 |t|,. ^-^ (p+»x) 2 , 
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^)] = E e * ^--jr) ( A7 ) 



where the correction to the free energy 6F is 

(p+4>x) 2 



V^e-^n A sin 2 (.Tr) [ 



l e -2«3:r(p+</,x)-n 2 T-T 2 -n 2 f rj ^g^ 



where actually j = 1/2MR 2 . At small r there are J cLt/t integrals and therefore a cutoff 1/oj c is needed. At low 
temperatures i < 1 one can retain only p = p' = and then the cutoff is not needed, as found below. Hence for 
t < 1, 

r/3/2 ^2-7^2 2 

6F = a S T a n —T7 -[l-e- n2 ^ r2 -™ 2 ¥ T cosh(2nr^T/t)] + 0(e^ 1/ *lnw c T) (A9) 

^ Jo sin (ttTt) 

The effective mass M* is defined in terms of the curvature, so that the 1st order correction is 

I d 2 8F ^ f 0/2 TT 2 T 2 2T ! 2 2T 

^ = w'° — ^ (2 " TT/i)V "^ (A10) 

Note that there is no divergence at r = 0. The dominant integration range is r < t/Tn 2 so that the 1st term in the 
exponent can be expanded; keeping terms to order t 2 we obtain in terms of x = rn 2 /2MR 2 , 

r M n ^ f°°. ir 2 t 2 2 t 2 t 2 4 

n u 

= -2a ]>>„(! - 4 + R- + 12)^1 + •••) (AH) 



Hence to 1st order in t 

^ = l-2a$>„+4 to £^ (A12) 

n n 

At t = this result is consistent with Eq. 9 of Ref. [l3l 

The following sum rules are useful for evaluating these sums. Integrating Eq. ()A1 1) J™ dz we obtain: 

f> = 2--f , ^ (A13) 

n=l / 4r 2 sin 2 i z + j 



Fourier transform of Eq. (|A1[) 



1 ^^^=i^^^= ] cosnzdz (A14) 

/4r 2 sin 2 |z + 1 



and performing the n summation, we obtain 



DC: 



E"„ 1 /' 1 , i ~ / 1 . 

^ = ZL < T~^ + T) dz - ( A15 ) 



o A /4r 2 sin 2 \z+ 1 
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